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Fig. 10.26 Circular aperture geometry. 



still smaller, e.g. the four corner peaks (whose coordinates 
correspond to appropriate combinations of $' = ±3tt/2 and 
a' = ±3ti/2) nearest to the central maximum, each have 
relative irradiances of (^y) J . 

10.2.5 The Circular Aperture 

Fraunhofer diffraction at a circular aperture is an effect of 
very great practical significance in the study of optical in- 
strumentation. Envision a typical arrangement; plane waves 
impinging on a screen I containing a circular aperture and 
the consequent far-field diffraction pattern spread across 
a distant observing screen a. By using a focusing lens L 2 , 
a can be brought in close to the aperture without changing 
the pattern. Now, if L 2 is positioned within and exactly fills 
the diffracting opening in I, the form of the pattern is 
essentially unaltered. The light wave reaching £ is cropped 
so that only a circular segment propagates through L 2 to 
form an image in the focal plane. But quite obviously this is 
the same process that takes place in the eye, a telescope, 
microscope or camera lens. The image of a distant point 
source as formed by a perfectly aberration -free converging 



lens, is never a point, but rather some sort of diffraction 
pattern. We are essentially collecting only a fraction of the 
incident wavefront and cannot therefore hope to form a 
perfect image. As shown in the last section, the expression 
for the optical disturbance at P, arising from an arbitrary 
aperture in the far-field case, is 

c 0 n<ot- km r r 
E . A ^ e wr9 + z# ! * dS. [10.41 ] 

Aperture 

For a circular opening, symmetry would suggest introducing 
spherical polar coordinates in both the plane of the aperture 
and the plane of observation, as shown in Fig. 10.26. 
Therefore, let 

z - p cos <p y = p sin <p 
Z = q cos <I> K = <7sin<I> 

and so the differential element of area is now 
dS = pdpd<p. 
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Substituting these expressions into Eq. (10.41), it becomes 
E = e /<fc M .R>coM*-e>,, dp d(t> (1 o 46) 

n Jp = 0 J«b = 0 

Because of the complete axial symmetry, the solution must be 
independent of <D. We might just as well solve Eq. (10.46) 
with <D = 0 as with any other value, thereby simplifying 
things slightly. 

The portion of the double integral associated with the 
variable <p. 

is one that arises quite frequently in the mathematics of 
physics. It is a unique function in that it cannot be reduced 
to any of the more common forms, such as the various hyper- 
bolic, exponential or trigonometric functions and indeed, 
with the exception of these, it is perhaps the most often 
encountered. The quantity 



J 0 (u) = 




(10.47) 



is known as the Bessel function (of the first kind) of order 
zero. More generally 

J m (u) = e i(«.^-co.r)^ (10.48) 

2n Jo 



represents the Bessel function of order m. Numerical values 
of J 0 (u) and J x (u) are tabulated for a large range of u in 
most mathematical handbooks. Just like sine and cosine, 
the Bessel functions have series expansions and are cer- 
tainly no more esoteric than these familiar childhood 
acquaintances. As seen in Fig. 10.27, J 0 (u) and J x (u) are 
slowly decreasing oscillatory functions that do nothing 
particularly dramatic. 

Equation (10.46) can be rewritten as 

f = 2n J J 0 (kpq/R)p dp. (10.49) 

Another general property of Bessel functions, referred to 
as a recurrence relation, is 

When m « 1, this clearly leads to 

J uUaii/) du 1 = uJ x (u) (10.50) 

with u' just serving as a dummy variable. If we now return 
to the integral in Eq. (10.49) and change the variable such 
that w = kpq/R, then dp = (R/kq) dw and 

I Jo{kpq/R)pdp=(Rfkq) 2 [ J 0 (w)wdw. 

Jp-0 Jh- = 0 
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Making use of Eq. (10.50), we get 
Z A e lic *~ kR> 



£(t) = 



2na 2 (R/kaq)J i (kaq/R). (10.51) 



The irradiance at point P is <(Re £) 2 > or \EE \ that is 
<? 2 L kaq/R J ' 



/ = 



(10.52) 



where A is the area of the circular opening. To find the ir- 
radiance at the center of the pattern, i.e. at P 0 , set q = 0. It 
follows from the above recurrence relation (m - 1 ) that 



J 0 (u) = -j-JtM + • 

du u 



(10.53) 

From Eq. (10.47) we see that J 0 (0) = 1 and from Eq. (10.48), 
Ji{0) = 0. The ratio of J x (u)/u as u approaches zero has 
the same limit (L'Hopital s rule) as the ratio of the separate 
derivatives of its numerator and denominator, namely 
dJ^/du over one. But that means that the right-hand 
side of Eq. (10.53) is twice that limiting value, so that 
J x (u)/u = \ at u = 0. The irradiance at P 0 is therefore 



/(0) = 



ZjA 2 
2R 2 



(10.54) 



1(0) 



which is the same result obtained for the rectangular opening 
(10.43). If R can be presumed essentially constant over the 
pattern, we can write 

Since sin 0 = q/R, the irradiance can be written as a func- 
tion of 0; 

pj(* a sin Q )T {1056) 

and as such, is plotted in Fig. 10.28. Because of the axial 
symmetry, the towering central maximum corresponds to a 
high-irradiance circular spot known as the Airy disk. It was 
Sir George Biddell Airy (1801-92) Astronomer Royal of 
England, who first derived Eq. (10.56). The central disk is 
surrounded by a dark ring which corresponds to the first 
zero of the function J^u). From the standard tables J x (u) 
= 0 when u = 3.83 i.e. kaq/R = 3-83. The radius drawn 
to the center of this first dark ring can be thought of as the 
extent of the Airy disk. It is given by 

q ^22~. (10.57) 
* l 2a 




kas'inO 



Fig. 10.28 The Airy pattern. 

For a lens focused on the screen a, the focal length / ^ R 
and so 



7,-1.22-, 



(10.58) 



where D is the aperture diameter i.e. D = 2a. (The diameter 
of the Airy disk in the visible spectrum is very roughly equal 
to the //# of the lens in millionths of a meter.) As shown 
in Figs. 10.29 to 10.31, q x varies inversely with the hole's 
diameter. As D approaches k, the Airy disk can be very large 
indeed and the circular aperture begins to resemble a point 
source of spherical waves. 

The higher-order zeroes occur at values of kaq/R equal 
to 7.02, 1 0.1 7, etc. The secondary maxima are located where 
u satisfies the condition 



d_ 
du 



which is equivalent to J 2 (u) = 0. From the tables then, 
these secondary peaks occur when kaq/R equals 5.14, 8.42, 
11.6, etc.; whereupon 1/1(0) drops from one to 0.0175, 
0.0042 arid 0.001 6 respectively. 



10.2 




Fig. 10.29 Airy rings (0.5 mm hole diameter). [Photo by E.H.] 




Fig. 10.30 Airy rings (1 0 mm hole diameter). [Photo by E H ] 

Circular apertures are preferable to rectangular ones as 
far as ler.s shapes go, since the circle's irradiance curve is 
broader around the central peak and drops off more rapidly 
thereafter. Exactly what fraction of the total light energy 
incident on a is confined to within the various maxima is a 
question of interest, but one somewhat too involved to solve 
here.* On integrating the irradiance over a particular region 



• See Born and Wolf, Principles of Optics, p. 398 or the very fine 
elementary text by Towne. Wave Phenomena, p. 464. 
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(b) 

Fig. 10.31 (a) Airy rings — long exposure (1 .5 mm hole diameter), 
(b) Central Airy disc — short exposure with the same aperture. 
[Photos by E.H.] 

of the pattern, one finds that 84% of the light arrives within 
the Airy disk and 91% within the bounds of the second dark 
ring. 

10.2.6 Resolution of Imaging Systems 

Imagine that we have some sort of lens system which forms 
an image of an extended object. If the object is self-luminous 
it is likely that we can regard it as made up of an array of 
incoherent sources. On the other hand an object seen in re- 
flected light will surely display some phase correlation be- 
tween its various scattering points. When the point sources 



